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ABSTRACT 



If aspherical dust grains are immersed in an anisotropic radiation field, their tem- 
^ ■ perature depends on the cross-sections projected in the direction of the anisotropy. It 

^^ ■ was shown that the temperatm'e difference produces polarized thermal emission even 

V/^ . without alignment, if the observer looks at the grains from a direction different from the 

anisotropic radiation. When the dust grains are aligned, the anisotropy in the radiation 
Q^ I makes various effects on the polarization of the thermal emission, depending on the 

^^ I relative angle between the anisotropy and alignment directions. If the both directions 

are parallel, the anisotropy produces a steep increase in the polarization degree at short 
wavelengths. If they are perpendicular, the polarization reversal occurs at a wavelength 



^ ■ shorter than the emission peak. The effect of the anisotropic radiation will make a 

c/3 . change of more than a few % in the polarization degree for short wavelengths and the 

effect must be taken into account in the interpretation of the polarization in the thermal 

emission. The anisotropy in the radiation field produces a strong spectral dependence 

r> ' of the polarization degree and position angle, which is not seen under isotropic radia- 

C^ I tion. The dependence changes with the grain shape to a detectable level and thus it 

will provide a new tool to investigate the shape of dust grains. This paper presents 

examples of numerical calculations of the effects and demonstrates the importance of 

anisotropic radiation field on the polarized thermal emission. 

Subject headings: dust, extinction — infrared: ISM: continuum — polarization 



1. INTRODUCTION 

The polarization observation in the far-infrared region is a quite efficient means for the investi- 
gation of magnetic fields in dense regions and the alignment mechanism of dust grains (Hildebrand 
1988). A large number of polarization observations have been reported recently for various in- 
teresting objects in the far-infrared and the submillimeter regions (e.g., Hildebrand et al. 1995; 



Novak et al. 1997; Schleuning 1998; Dowell 1997; Dowell et al. 1998). Hildebrand et al. (1999) have 
provided a summary of the latest observations and discussed several factors which could affect the 
observed polarization and the spectral dependence, such as the different populations of grains with 
different temperatures, other than the magnetic field and grain alignment. In this paper, another 
mechanism, the anisotropy of the incident radiation field, which is believed to occur commonly in 
celestial objects, is shown also to affect the polarization of thermal emission. 

If aspherical grains are immersed in an anisotropic radiation field, the grains with larger cross- 
sections projected in the radiation direction absorb more photons and thus have higher temper- 
atures than those with smaller cross-sections. The former grains emit the radiation polarized in 
the direction perpendicular to the incident radiation field, while the radiation from the latter is 
polarized in the parallel direction, if one observes them from a direction different from the radia- 
tion anisotropy. Therefore owing to the difference in the temperature, the resultant total thermal 
emission is polarized perpendicularly to the incident radiation. 

Onaka (1995, hereafter Paper I) has made numerical calculations of this effect and shown 
that it produces polarization of a detectable degree (~ 1%) around 100 ^m and a steep rise in the 
polarization degree at shorter wavelengths even without any net alignment. This paper presents 
the numerical results of several simple examples when aligned dust grains are immersed in an 
unidirectional radiation field and discusses the effects of anisotropic radiation on the polarized 
thermal emission from aligned dust grains. In §2, the basic assumptions made in the present 
study are described, and the numerical results of several examples are presented in §3. In §4, the 
results are examined in relation to the observations and the effects of anisotropic radiation on the 
polarization of thermal emission are discussed. 



2. BASIC CONFIGURATION 

Figure 1 depicts the geometry adopted in this paper: The radiation field is assumed to be 
unidirectional and coming from -|-Z direction; The body axis of greatest inertia of the grain A is 
tilted from the Z axis by ip; The angle between the alignment axis Z' and the Z axis is denoted by 
7 and the Z' axis is assumed to be in the Y-Z plane; The observer is located in the X-Z plane with 
the angle /? from the X axis. In the following, /3 is set as zero (the observer is on the X axis) for 
simplicity. 

Dust grains perform complex motions in the interstellar space; the grain axis A makes a 
precession about the angular momentum J and J makes a precession about the magnetic field due 
to the magnetic moment of the grain (Lazarian 1997b). The former angular velocity is of the order 
of the grain angular velocity. The Barnett effect, the inverse of the well-known Einstein-de Haas 
effect, produces the magnetic moment of the rotating grain, which is much lager than that acquired 
by a charged rotating grain (Dolginov & Mytrophanov 1976). The period of the precession due to 
the Barnett effect magnetic moment is of order of 10^ s (Lazarian 1994). It is still much longer than 
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Fig. 1. — A schematic diagram of the assumed configuration. In the numerical calculations pre- 
sented in this paper /3 is set as zero. 



the mean duration time of photon absorption in the interstellar radiation field (ISRF) of 1 — 10 s 
for 0.01 — O.l/um-sized grains (Purcell 1976). Therefore, the angular momentum direction can be 
assumed to be fixed in the reference frame as long as the temperature of the grain is considered. 

Dust grains in the interstellar space are thought to rotate suprathermally due to the irreg- 
ularities in the shape or surface characteristics (Purcell 1979; Draine & Weingartner 1996). The 
rotation time scale is very short, typically 10~^ to 10~^ s, depending on the grain size, and on 
the mechanism and efficiency of suprathermal rotation under the physical conditions in question. 
Owing to the internal energy dissipations, such as the anelasticity or the Barnett effect relaxation, 
the rotation axis aligns with the body axis of greatest inertia if the grain rotates suprathermally 
(Purcell 1979). Most recently Lazarian & Draine (1999b) have shown that nuclear spin relaxation 
within the grain will be the most efficient process in the internal energy dissipations. On the other 
hand, Lazarian (1994) pointed out that the alignment is always not perfect and the deviation should 
be correctly taken into account when the internal temperature is comparable with the grain kinetic 
temperature (Lazarian 1997a; Lazarian & Roberge 1997). The deviation is, however, typically very 
small for suprathermally rotating grains (Lazarian 1994). In this paper, the body axis A is assumed 
to coincide with the angular momentum J for simplicity. 

The alignment mechanism of interstellar dust grains is still a hot issue. Various mechanisms 
and their revisions have so far been proposed, including paramagnetic relaxation of ordinary para- 
magnetic grains (Davis & Greenstein 1951) or grains with superparamagnetic inclusions (Jones & 



Spitzer 1967; Mathis 1986; Martin 1995; Goodman & Whittet 1995), mechanical alignment (Gold 
1951; Lazarian 1994), alignment due to suprathermal rotation (Purcell 1979; Lazarian 1995a, b; 
Lazarian & Draine 1997, 1999a), and alignment by radiation torques (Harwit 1970; Dolginov & 
Mytroplianov 1976; Draine & Weingartner 1996). The efficiency of each mechanism depends on the 
physical conditions in question. See Lazarian, Goodman, & Myers (1997) for a latest review and 
details of each mechanism. The time scale of paramagnetic relaxation and mechanical alignment 
is, however, typically much longer than the photon incident duration mentioned above. The mean 
duration of the crossover for suprathermally rotating grains must also be quite long if this mecha- 
nism is effective in the grain alignment. In this paper, I assume that the grain has a fixed rotation 
axis, which coincides with the body axis of greatest inertia, relative to the incident radiation and 
thus the temperature of the grain is given by a function of ip. 

The alignment distribution of the dust grains around the alignment axis depends on the mech- 
anism of alignment. In this paper, I assume the following distribution function for simplicity: 

fie)dcose = - — . 2a\-J2 dcosO, (1) 

(cos^ & + q sm t))^/'^ 

where g is a parameter of the alignment degree. When no alignment is established, q is equal to 
unity. For q > 1, the distribution peaks at ^ = 0, while the perpendicular alignment is expressed 
hy q < 1. In the following, only the cases of q > 1 are considered. 

For the Davis-Greenstein alignment, the distribution is given by equation (1) and q is equal 
to Tg/Tav, where Tg is the gas temperature and Tav is the weighted mean of gas and dust tem- 
peratures (Jones & Spitzer 1967). The distribution with q = Tg/Tav may be applicable only if 
grains are superparamagnetic or ferromagnetic. In the mechanical alignment by supersonic flow 
the distribution is given by the same functional form with q = (1 + g)/{l + g + s), where 5 is a 
parameter of the grain shape and s is that of the flow anisotropy (Lazarian 1994). It can also 
be shown that equation (1) holds for the alignment by suprathermal rotation with q being equal 
to exp(2tc/tr), where tc and tr are the time from the last crossover and the magnetic time scale, 
respectively, if the orientation after a crossover is completely isotropic (Purcell 1979). Lazarian & 
Draine (1997, 1999a) have shown that the angular momentum before and after a crossover should 
be correlated. This effect may be taken into account by a slight modification in the parameter 
(Lazarian & Draine 1997). When the coincidence of the body axis and the angular momentum 
axis is not perfect, further modification of the distribution is needed (Lazarian 1997a; Lazarian & 
Roberge 1997). The functional form of equation (1) may not always hold. It is adopted here only 
to represent a simple alignment distribution. 

According to Lazarian & Draine (1997), the grain alignment may be nearly perfect. In the 
perfect alignment case, the present results may still be applicable to the case that the alignment 
direction (the angle 7) changes along the line of sight and the distribution of 7 is given by a function 
similar to equation (1) provided that the direction of the radiation field is fixed. 



The measure of alignment a is defined by (Greenberg 1968) 
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where <> denotes the distribution average. Taking the ensemble average with equation (1), one 
obtains (Purcell 1979) 
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Figure 2 plots the distributions given by equation (1) with q = 1.2, 1.5, and 2.0. These 
correspond to a = 0.037, 0.083, and 0.14, respectively. It can be seen that the distribution is quite 
broad and the alignment is far from being perfect for this range of q . 
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Fig. 2. — The alignment distribution function f{9) given by equation (1). The dotted line indicates 
the distribution for q = 1.2, the dashed line for q = 1.5, and the solid line for q = 2.0. 

The shape of the interstellar grains is not well known. Several past investigations have sug- 
gested that oblate grains fit the observations better than prolate grains (Lee & Draine 1985; Henning 
& Stognienko 1993; Kim & Martin 1995; Hildebrand &: Dragovan 1995). Hildebrand & Dragovan 
(1995) have indicated that oblate spheroidal grains with axial ratio b/a = 1.5 fit the mid-infrared 
observations best, while Kim &: Martin (1995) have suggested that flat oblate grains fit the visual 
polarization observations most satisfactorily. Kim &: Martin (1995) have also indicated that dust 



grains smaller than 0.1 /im are not well aligned. On the other hand, alignment of silicate grains 
are indicated by the observed mid-infrared band feature in the polarization spectrum (Aitken et al. 
1988). 

In the following calculations, the grains are assumed to be oblate silicate spheroids. The 
calculations are made with the axial ratio h/a = 1.5,2.0, and 3 to investigate the dependence on 
the grain shape. The particle size is fixed as a = 0.1 /^m. The optical properties are taken from 
those of astronomical silicate (Draine & Lee 1984) and the incident radiation is adopted from the 
ISRF given by Mathis, Mezger, & Panagia (1983). The far- infrared emission is assumed to be 
optically thin. The absorption cross-sections of spheroidal grains are calculated by the method of 
expansion series of spheroidal functions developed by Onaka (1980). Details of the formulation 
have been given in Paper I. 



NUMERICAL RESULTS 
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Fig. 3. — The temperature of oblate silicate spheroid grains of a = 0.1 ;um in a unidirectional 
radiation field against the cosine of tilt angle '(/'• The solid line indicates the grain with hja = 1.5, 
the dashed line that with hja = 2, and the dot-dashed line that with hja = 3. 



The temperature of the grain is plotted against cos'^ in Figure 3 for the parameters described 
in the previous section. The temperature of the grain with small axial ratio is higher than that 
with large ratio because the size of the major axis a is fixed in the present calculation and thus the 



total volume is small for the grains with small axial ratio. The difference in temperature between 
cos'i/'=0 and 1 is quite small even for the oblate spheroid of axial ratio of 3. It becomes smaller 
as the grain becomes less flattened. The effect is insignificant in the Rayleigh-Jeans tail of the 
thermal emission, but becomes appreciable for shorter wavelengths than the emission peak (see 
Paper I). In this section, three cases of the direction of alignment relative to the radiation direction 
are investigated. 
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Fig. 4. — a) Spectral dependence of the polarization degree (left axis) for the alignment axis parallel 
with the radiation direction (7 = 0°) for the oblate spheroidal grain of b/a = 2 together with the 
emission spectrum (right axis in arbitrary units; dotted line labeled as flux). The dashed lines 
indicate the polarization degrees expected for the isotropic radiation field, while the solid lines are 
those for a unidirectional radiation field. The numbers show the values of alignment parameter q. 
b) The same as a) but for different grain shapes. The solid lines indicate the results for b/a = 3, 
the dashed lines those for b/a = 2.0, and the dotted lines for b/a = 1.5. The thick lines are for 
q = 1.5 and the thin lines are for q = 2.0. 

Figure 4a shows the spectral dependence of the polarization degree (left axis; solid lines) of 
the oblate grain of b/a = 2 for 7 = 0° together with the emission spectrum (right axis; dotted line 
labeled as flux). The results of isotropic radiation case are plotted together by dashed lines for 
comparison. As shown in Paper I, even without alignment {q = 1), thermal emission of the dust 
grains in an anisotropic radiation field is polarized and the polarization degree increases steeply as 
the wavelength becomes shorter. For the cases when the grains are aligned, the anisotropic radiation 
field increases the polarization degree appreciably at short wavelengths compared to the isotropic 
radiation case. Under the isotropic radiation field, almost no spectral dependence is expected for 
wavelengths longer than 30 ^m (Hildebrand 1988). The strong spectral dependence due to the 



anisotropic radiation field is quite remarkable. The effect of anisotropic radiation field becomes 
more than a few % in the polarization degree particularly for wavelengths shorter than the thermal 
emission peak. This can easily be understood from the nature of the mechanism. 
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Fig. 5. — Spectral dependence of the polarization degree for the alingment axis perpendicular to 
the radiation direction (7 = 90°). The solid lines indicate the grain with b/a = 2.0; the thin line 
for q = 1.2, the medium-thickness line for q = 1.5, and the thick line for q = 2.0. The dashed 
line shows the result for the grain of b/a = 1.5 and q = 1.5, and the dotted line for b/a = 3.0 and 
q= 1.5. 

The spectral dependence of the polarization degree changes with the grain shape (Figure 4b). 
Even when the polarization degree is similar at longer wavelengths (e.g., the grain of b/a = 1.5 
with q = 2.0 (thin dotted line) and that of b/a = 3 with q = 1.5 (thick solid line)) the flatter 
grain shows a steeper spectral dependence. This is because the effects of anisotropic radiation are 
larger for flatter grains than less flatter grains due to the larger difference in the cross-sections. Less 
alignment degree is required for flatter grains to produce a given polarization degree at wavelengths 
longer than 100 fim. 

When the alignment axis is perpendicular to the incident radiation direction (7 = 90°), the 
polarization degree of aligned grains first decreases as the wavelength becomes shorter (Figure 5). 
Then the polarization reversal occurs at some wavelength where the effect of anisotropic radiation 
starts dominating over the original alignment. In the calculation, the sign of the polarization degree 
is defined as positive if the emission is polarized in the Y axis. The reversal wavelength depends on 



the alignment degree and grain shape, and occurs at a wavelength shorter than the emission peak. 
Figure 5 also indicates a similar trend of the spectral dependence with the grain shape to that for 
7 = 0°. 
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Fig. 6. — Spectral dependence of the polarization degree for 7 
as in Figure 5. 
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If the alignment axis is tilted against the radiation direction by an arbitrary angle, further 
complicated effects appear. Figures 6 and 7 plot the results for 7 = 45°. The spectral dependence 
of the polarization degree is similar to that for 7 = 0° or 90°, while the polarization angle also varies 
with the wavelength (Figure 7). The polarization angle is parallel with the alignment direction for 
long wavelengths (45°), but it becomes perpendicular to the radiation direction as the wavelength 
becomes shorter. The spectral dependence of the polarization angle also varies with the grain shape 
and alignment degree. For low alignment degree case (g = 1.2, thin solid line) the deviation from 
the original alignment direction can be seen even for wavelengths longer than 100 fira. 



DISCUSSION 



Hildebrand et al. (1999) have presented the latest results of far- infrared polarization measure- 
ments and shown that most of the observed polarization degrees in 60 — 100 fiui are in the range 
less than 10%. The anisotropic radiation field produces very small differences in temperature of 
the dust grains (Figure 3), but the effect in the polarization degree is still a few percent level for 
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Fig. 7. — Spectral dependence of the polarization angle for 7 
in Figure 5. 
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wavelengths shorter than 100 /xm, being compatible with the observed range. Thus the effect of 
anisotropic radiation field cannot be neglected in the interpretation of the polarization observation 
of thermal emission. 

Hildebrand et al. (1999) have also shown the spectral dependence of the polarization degree in 
the far-infrared region. In the full samples presented, there is a general trend in the distribution of 
observed degrees of polarization that the degree decreases with the wavelength, though this trend 
is seen in the combined data set obtained at different regions and/or objects. In two particular 
objects the polarization degree either increases as the wavelength decreases (Orion core) or decreases 
sharply for short wavelengths (M17). In another case Dowell (1997) showed that the polarization 
degree increases from submillimeters to 60 /xm in Sgr B2 core and claimed that the polarization 
in the far-infrared in Sgr B2 core is due to selective absorption. In most other observed objects, 
however, the far-infrared optical depth is not very large and the polarization occurs in emission. 
Hildebrand et al. (1999) attributed the observed spectral dependence to two populations of grains 
that differ in their polarization efficiencies and in their temperatures due to the difference in their 
optical properties. 

The mechanism proposed in the present paper also comes from the difference in the tempera- 
ture of dust grains, but it originates from the anisotropy in the incident radiation field, which is not 
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uncommon in most astronomical objects. It could make either positive or negative spectral depen- 
dence as being actually observed, depending on the configuration of the alignment and radiation 
field. The change in the polarization angle could also occur due to this mechanism. In the FIR 
polarization observations of Sgr B2, the observed change in the polarization angle from the core to 
envelope was successfully explained by the selective absorption in the core region, where the optical 
depth is expected to be larger than unity (Dowell 1997; Novak et al. 1997). The present mechanism, 
on the other hand, can produce polarization angle variation even in optically-thin regions. It should 
be noted that in the vicinity of the bright source the anisotropic radiation could also be effective 
in the grain alignment (Draine & Weingartner 1996). 

Because of the nature of the blackbody radiation, the effect due to the difference in temper- 
ature appears most severely at short wavelengths. The present study has in fact shown that the 
polarization at long wavelengths is not affected largely and thus indicates the alignment character- 
istics of the dust grains directly. The information such as the direction of the magnetic field may 
be most correctly traced by polarization observations for wavelengths longer than 100 fiui, though 
spectral observations should be quite important to investigate the effects other than the alignment. 

On the other hand the polarization at short wavelengths should contain the information of the 
dust properties. It is generally difficult to separate the effects of the grain shape and alignment 
degree solely from the polarization observation of continuum emission. However, the effect of 
radiation anisotropy may provide a possible method to separate these factors. If we know the 
radiation geometry (location of the heating source) quite well by other observations, the effect can 
be estimated to some accuracy. As shown in Figures 4, 5, and 6, the polarization degree differs 
at short wavelengths for grains with different shapes even if the degree at long wavelengths is the 
same. A larger polarization degree is expected for flatter oblate grains at short wavelengths and the 
differences in the degree among the grains of b/a = 1.5,2, and 3 will be in a distinguishable level 
(~ 1%). If the alignment and radiation directions are not parallel, the change of the polarization 
angle should provide an independent confirmation of the derived parameters. Hence, the effect of 
anisotropic radiation field may provide a new tool to investigate the grain shape other than the 
dust band (Lee & Draine 1985; Aitken et al. 1988; Henning & Stognienko 1993; Hildebrand et al. 
1995), though in an actual application one may also have to take account of the effects of dust 
populations and/or the optical depth (Hildebrand et al. 1999). 

The effect of particle size may also have to be examined. It is coupled with the spectrum of the 
incident radiation field (Paper I) and has to be investigated with a radiation spectrum appropriate 
for the object in question. For the ISRF spectrum, the effect of particle size is not significant between 
a = 0.1 /um and 0.03 nm because its energy peaks at a sufficiently long wavelength. However, in the 
case where high energy photon dominates, smaller grains can produce larger polarization (Paper 

I)- 

The present study demonstrates the importance of the spectral observations of polarization 
in the far-infrared. Correct information on the alignment or magnetic field may require careful 
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polarization observations over a range of wavelengths. Spectral observations of the polarization 
degree and angle in the thermal emission should be crucial for the investigation of the effect of the 
anisotropic radiation field. The present effect may also be useful in the identification of the heating 
source for very obscured objects. 

This work is stimulated by the ISO polarization workshop held at Vilspa ISO Data Center. 
The author would like to thank the workshop organizers for their efforts in the arrangement of the 
ISO polarization workshop. This work is supported by Grant-in-Aids for Scientific Research from 
Japan Society for the Promotion of Science (No. 10559005 and 11440062). 
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